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Abstract 

In a model of 3-brane embedded in 5D space-time we calculate the graviton emission from the 
brane to the bulk. Matter is confined to the brane, gravitons produced in reactions of matter 
on the brane escape to the bulk. The Einstein equations which are modified by the terms 
due to graviton production are solved perturbatively, the leading order being that without the 
graviton production. In the period of late cosmology, in which in the generalized Friedmann 
equation the term linear in the energy density of matter in dominant, we calculate the spectrum 
of gravitons (of the tower of Kaluza-Klein states) and the collision integral in the Boltzmann 
equation. We find the energy-momentum tensor of the emitted gravitons in the bulk, and using 
it show that corrections due to graviton production to the leading-order terms in the Einstein 
equations are small, and the perturbative approach is justified. We calculate the difference of 
abundances of 4 He produced in primordial nucleosynthesis in the models with and without the 
graviton production, and find that the difference is a very small number, much smaller than 
that estimated previously. 



1 Introduction 



Brane-world scenarios with a 3-brane identified with the observable Universe which is embedded in 
a higher-dimensional space-time provide an alternative to the standard 4D cosmology [1, 2, 3, 4, 
5], reviews [6, 7, 8]. A necessary requirement on these models is that they should reproduce the 
main observational cosmological data, the age of the Universe, abundances of elements produced 
in primordial nucleosynthesis, etc. A general property of the models with extra dimensions is that 
gravity propagates in the extra dimensions independent of whether the ordinary matter is confined to 
the brane or not. This entails a peculiar property of the models with extra dimensions which is absent 
in the standard cosmology: gravitons which are produced in reactions of particles of matter on the 
brane can escape from the brane and propagate in the bulk [9, 10, 11, 12, 14, 15]. As a consequence 
the Einstein equations contain terms accounting for the graviton emission. Cosmological evolution 
of matter on the brane is also affected by this process. 

Roughly the energy loss on the brane due to the process a + b — > G + X can be estimated as 

-tt = - < n a n b a a+b ^ G+x vE G >, 
at 

where in the radiation-dominated period of the evolution of the Universe n a ,nb ~ T 3 and Eq ~ T. 
This yields 

dt 
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Here T is temperature of the Universe and k 2 = 8ir/M 3 , where M is the 5D Planck mass, is the 5D 
gravitational constant. Although at the late stages of the cosmo logical evolution the energy loss to 
the bulk is small, at high temperatures this effect can be significant. 

In the present paper we consider the problem of graviton production and emission to the bulk 
in a model with one 3-brane embedded in the bulk with one extra dimension. We perform our 
calculations in a picture in which the metric is time-dependent and the brane is located at a fixed 
position in the extra dimension. 

The Einstein equations including the terms due to graviton production are solved in the pertur- 
bative approach. In the leading order we neglect the graviton production, and include it in the next 
order. The classical background metric is the leading-order solution of the 5D Einstein equations. 
Background metric is a warped extension of the metric on the brane to the bulk with the warping 
factor e ±m [2]. Here \i ~ y/—A, where A/k 2 is the 5D cosmological constant and y a coordinate of 
the extra dimension. 

We perform calculations in the period of late cosmology. In this period in the generalized Fried- 
mann equation 



the term linear in energy density on the brane is much larger than the quadratic term, or, equivalently, 
/i ^> K 2 p. For the value of \i of order 10~ 12 GeV, which we adopt for numerical estimates, the range of 
temperatures of the Universe characteristic to the period of late cosmology extends to T ~ 5-10 2 GeV. 

Evolution of the energy density on the brane p is determined by the Boltzmann equation. The 
collision term in the Boltzmann equation accounts for the graviton emission from the brane. To 
calculate the collision term, solving the equations of motion for fluctuations over the background 
metric, we find the spectrum and the eigenfunctions of of the tower of Kaluza-Klein gravitons and 
the energy-momentum tensor of the emitted gravitons. 

Solution of the Einstein equations requires integration of the 5D energy-momentum tensor over 
coordinates of the extra dimension. This, in turn, requires determination of the form of the 5D 
energy-momentum tensor not only on the brane, where gravitons are emitted, but also in the bulk. 
Integrating the geodesic equations for the null geodesies along which gravitons propagate and using 
conservation equations for the energy-momentum tensor in the bulk, we find the form of the energy- 
momentum tensor in the bulk. 

Using the obtained energy-momentum tensor in the bulk, we show that in the period of late 
cosmology in the Friedmann equation corrections to the leading-order terms due to graviton emission 
are small, and the perturbative approach is consistent. 

Graviton emission changes cosmological evolution of matter on the brane. Time (temperature) 
dependence of of the Hubble function determined from the Friedmann equation which includes the 
terms due to graviton emission is different from that in the standard cosmological model. This, in 
turn, results in a change of abundances of light elements produced in primordial nucleosynthesis [16]. 
We find that the difference of abundances of 4 He produced in primordial nucleosynthesis calculated 
in the models with and without the graviton production is a small number, much smaller than that 
estimated in [10, 11, 12]. A crude estimate of the graviton production in the early cosmological 
period, which takes ito account the bounce of produced gravitons back to the brane [12], does not 
alter this result. 




9 



2 System of Einstein equations 

We consider the 5D models with one 3D brane embedded in the bulk with the action 



1 

2^ 



J d b x^-g^(R^ - 2A) + 2 j K - j d 4 x^-g^) a- J d 4 x^-g^L m , (1) 
.s as J 9S dT, 



where x 4 = y is coordinate of the infinite extra dimension, k 2 = 8n/M 3 . We consider a class of 
metrics of the form 1 



ds'l = gffdx'dx 3 = -n' 2 (y,t)dt 2 + a 2 (y,t)n ab dx a dx b + dy 2 = dy 2 + g^dx^dx" . 



(2) 



The brane is spatially flat and is located at the fixed position y — 0. The 5D Einstein equations are 

J^g~W 

9ij^ ~ C / — 5 (y)9»vK a. (3) 



dj — Rij ^9ijR — k 2 T^ ^ 



Here is the sum of the energy-momentum tensor of matter confined to the brane and the bulk 



energy-momentum tensor. The energy-momentum tensor of matter on the brane is taken in the form 



<%) t m = dia 9^(y){-p,P,p,p}- 

The components of the Einstein tensor G^- 
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satisfy the relations 2 (cf. [2]) 



n a a 
n a a 



G - G 4 -^ 
G44 — G04 — 



(4) 

(5) 
(6) 
(7) 



2a'a 3 
3 

2da 3 



F' 



F, 



where 



F = (da 



(aa) 
n 2 



(8) 
(9) 

(10) 



In the leading order we neglect emission of the gravitons produced in collisions of particles of matter 
on the brane in the bulk. In the next order we take into account the graviton emission into the bulk, 
the energy- momentum tensor in the bulk is Tj. 

1 The indices i, j run over 0,...,4, the Greek indices are 0,...,3, and a, b = 1, 2, 3 
2 Prime and dot denote differentiation over y and t. 
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For the following it is convenient to introduce the normalized expressions for energy density, 
pressure and cosmological constant on the brane which all have the same dimensionality [GeV] 



A n 2 a K 2 b K 2 p 

The functions a(y,t) and n(y,t) satisfy the junction conditions on the brane [2] 3 

a'(0,t) 



a(0,t) 
n'(0,t) 
n(0,t) 



= -a-p(t) (12) 
= 2p(t) + 3p(t) - a 



Using reparametrization of t, we set n(0,t) = 1, i.e. t is the proper time on the brane. Eq. (9) can 
be rewritten as 

F = - M V) - y (« 4 )'T 04 + y (a 4 )T 44 . (13) 



On the brane, using junction conditions and setting a ~ p, we have 
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F = /iV) + -^-(p + /i)T 04 + — 3—^44. (14) 

Integrating (14) in the interval (t, ti), where the initial time ti (of the onset of the period of the late 
cosmology) is defined below, we obtain 

t t 
F(0,t) = p 2 a\0,t) + J dt'f 04 (t')(p(t') + p)a 4 (0,t') + J dt'f 44 (t')a(0,t')a 3 (0,t') - C. (15) 

ti t t 

Substituting expression (10) for F and using the junction conditions, we rewrite (15) in a form 
of the generalized Friedmann equation (cf. [1, 2]) 

t 

H\t) = p\t) + 2pp(t) + pp w (t) - J£jL-^ J dt> [f 04 (t')(p(t') +p)+ f u (t')H(t')} a 4 (0, t>). (16) 

U 

Here we introduced 

d(0,t) , , /a(0,t ) x 1 



H ^ = Mry pM= "^{Wf) (17) 

where p w (t) is the so-called Weyl radiation term [7], pp w o = C/a A (0,t ). To comply with the 
observational data, the Weyl radiation term, which appears in the leading-order solution to the 
Einstein equations should be small 5 . Assuming that the terms with the bulk energy- momentum 



3 We assume invariance y -H- —y. 

4 From the 5D Einstein equations and the fit of the data on abundance of 4 He produced in nucleosynthesis it follows 
that cr 2 = [i 2 {l + 0(Hq/ ' ^ 2 )) 1 where -ffo is the present-time Hubble parameter [18, 19]. For \i ~ 10~ 12 GeV correction 
is - 1CT 60 . 

5 The Weyl radiation term p w is small as compared to the radiation energy density p r . BBN constraints to the ratio 
of the Weyl radiation term to the photon energy density at the temperatures T < 0.8MeV are —0.4 < p w /p~, < 0.1 
[22], or \ Pw / Pl \ < 0.55 [19]. 
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tensor can be treated perturbatively, we neglect corrections to the Weyl radiation term due to 
graviton emission from the brane. 

Substituting from the junction conditions the expressions for n'/n and a' /a at y = 0, we transform 
the (04) component of the Einstein equations (7) to the form 

p + 3H(p + p) = ^. (18) 

On the other hand, the same equation, which is a generalization of the conservation equation for the 
energy-momentum tensor of the matter confined to the brane to the case with the energy-momentum 
flow in the bulk, is obtained by integration across the brane of the 5D conservation law V«Tq = [9]. 

Some relevant equations following from general approach to the 3-brane embedded in the AdS 
bulk are found in the appendix A. 



3 Period of late cosmology 



Graviton production by hot matter is sufficiently intensive in the radiation-dominated period of 
cosmology. In the period of late cosmology the terms linear in energy density are dominant in the 
Friedmann equation. Equivalently, in the radiation-dominated period of late cosmology, this can be 
stated as p r (t)/p < 1, or 



p r (T) _ K 2 p r (T) _ 47r 3 ^(T)T 4 ^ 47r 3 ^(T)T 4 ^ i 



p 



6p 



90/iM 3 



90(pM pl y 



(19) 



where we used that pM 2 jM 3 ~ 1 [18, 19]. Taking p ~ 10 12 GeV, we find that the approximation 
of late cosmology is valid up to the temperatures of order 5 • 10 2 GeV. 

In the leading approximation, in which the bulk energy-momentum tensor is set to zero, from 
the system of the 5D Einstein equations and junction conditions on the brane is obtained the unique 
extension of the metric the components from the brane to the bulk [2] 



a 2 (y,t) = 



a 2 (0,t) 




(20) 



P P +2 + ^ 
p \p J p 



and 



n(y,t) 



a(0,t) 



Without the Weyl radiation term the function a 2 (y,t) is 



a 2 (y,t) 



a 2 {0,t) 



Mv\ ( E. 
p 



+ e 



-Mv\ ( E_ 
p 



+ 2 -2 



PIP 



p \p 



+ 2 



(21) 



(22) 



where a(0,t) is a solution of the Friedmann equation on the brane. a 2 (y,t) (22) As a function of 
y a 2 (y,t) (22) has the minimum equal to zero at the point \y\ 



P 



(23) 



°The expression for a 2 (y,t) in [2] contained no brane tension a . 



The function n(y, t) has no zeroes. At the point y the scalar curvature is finite. The function 
(20) has no zeroes, but the function n(y, t) has a zero. Again this zero is a coordinate singularity. 
In the region < \y\ < \y\ and for < 1 the functions a(y,t) and n(y,t) can be approximated 

as 

a(y,t) ~ a(0,t)e~ M (24) 
n(y,t) ~ e-"!"'. 

With (24) we obtain the approximate form of the metric in the region < \y\ < \y\ 

ds 2 ~ dy 2 + e~ 2M {-dt 2 + a 2 (0,t)r] ab dx a dx b ). (25) 

With the approximate metric we have a'(0, t)/a(0, t) = —p, which for p//i <C 1 is close to the exact 
junction condition a'(0,t)/a(0,t) = — (p + p). 

In the radiation-dominated period solution of the Friedmann equation 

H 2 = 2pp + p 2 (26) 

for the radiation energy density is 



1 

8pt r + At 

In the period of late cosmology the leading term is 

1 



p(t) = ^r—r,- ( 27 ) 



P(t)- 



8pt 2 

If p(t) > p, we have p(t) ~ 1/4*. 



4 Fluctuations of background metric 



To calculate the spectrum of the Kaluza-Klein tower of gravitons, we solve the field equations for 
fluctuations hy over the background metric g^. The purpose of this section is to clarify the issue of 
the gauge conditions in solving the equations for fluctuations. We show that the traceless transverse 
conditions on hij cannot be imposed as gauge conditions at the level of the action, but can be imposed 
on solutions of the equations of motion. 

The part of the action quadratic in fluctuations is 



,(5) 



(fl-A (5) ) 



R\h)h + 2R\h)K k 



(28) 



+- {2h qi , k h^ - h ik . q h ik ><> + h q h« - 2hjh ik ik ) 



The action I is invariant under the gauge transformations 

Ki = h kl - (V fc 6 + Vj^fc), 

where V is defined with respect to the metric g^. With the background metric (2), ds 2 
g^ u dx^dx u , the gauge transformations (29) take the form 

/144 = /144 — 20^C,4 



(D^ U + D U ^-2TU 4 ), 



(29) 
dy 2 + 

(30) 
(31) 
(32) 



where is defined with respect to the metric g^ u . 

To solve the equations following from the action (28) in the region < y < y, instead of the exact 
metric, we use the approximate metric (25). In the region < y < y we rewrite it as 

ds 2 = dy 2 + b(y)g lu ,(x)dx' t dx v , (33) 

where b(y) = e~ 2 ^. We impose the gauge conditions h Ai = 0. The gauge condition /i 44 = can be 
realized by performing transformation with 

Uv^) = \f dy'h u (y',x») (34) 

Condition = is imposed by taking 

Uy^ p ) = b(y) (c^x?) + £ dy'b-\y')(h 4 , - d^ 4 )(y',x^ . (35) 

Cfj,(x) are arbitrary functions. 

From the Einstein equations for the background metric (3) follow the relations 

1 D 2 

R = -A - -K 2 5(y)(r» - ASi) ~ 20/i 2 + I6p5(y) (36) 

R; = l -5^R - 2A) + K 2 5(y)(r; - ~ ^(-V + 2p5{y)) + K 2 5(y)r;, (37) 

where in the radiation-dominated period we have set = 0. The Ricci tensor calculated with the 
metric (33) is 

R; = -f (j + Q + R;(g) = ^(-V + 2n5{y)) + R;(g), (38) 

where we used the relations b'/b = — 2/isign(?/) and b"/b = 4/i 2 — Ap,5(y). Prime is derivative over y. 
Using (27), we obtain the estimate of the components of the tensor R^(g) 

r;(- 9 ) ~ i/t 2 ~ ^. 

The components of the tensor k 2 t^ are of order p. Comparing the expressions for R v calculated with 
the exact and approximate metrics, we find that in the period of late cosmology when p/p < 1 and 
/i 2 > Rp(g), p > k 2 t^, both expressions for the Ricci tensor coincide up to the terms of order p/p, 
which is the precision of our calculations in the period of late cosmology. 

Written with the approximate metric, the equations of motion are considerably simpler than with 
the exact metric. The (44) and (4p) and contracted {pv) components of equations of motion which 
follow from the action (28) are (cf. [17]) 

3 ft' 

V - D 2 h - —h! = 0, (39) 
(D„h - D'h,,,,)' = (40) 
D»D»h,„ - —ti - h" = 0, (41) 

where h = h flu g> J - u . From the Eqs. (39) and (41) it follows that 

b' 

h" + -h! = 0, 
o 
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or h'(y,x) = C(x)b~ 1 (y), where C(x) is an arbitrary function. Because of the reflection symmetry 
x -H- —x, the metric components g^ u (y,x) and h^ u (y,x) are even functions of y, and thus h'(y,x) is 
an odd function of y. Because b(y) is even in y, it follows that C(x) = and h'(y, x) = 0. We obtain 
that on the equations of motion h(y,x) is independent of y, h(y,x) = h(x). 

After imposing the gauge conditions = (35), there remain residual gauge transforma- 
tions ^(y,x) = b(y)C IJi (x). Under these transformations the trace h(x) transforms as 5h(x) = 
g^ u (D^C u (x) + D u C^(x)), and by a suitable choice of C^(x) can be transformed to zero. 

From the Eq. (40) with h — it follows that D^h' (y,x) = 0, or D /J, h lxl/ (y,x) = k u (x). Using 
the remaining residual gauge transformations with the functions C^(x) such that D^C^ = 0, we can 
impose the condition D tJ, h IJil/ (y,x) = 0. 



5 Eigenmodes and eigenvalues 

In the background of the approximate metric (33), in the gauge D^h^ u (y,x) = 0, M(y,x) = 0, the 
(fj,u) components of the field equations for fluctuations are 

h% - V V + b-\y)D p DPh^ + %)4/iV = 0. (42) 

Expanding the functions h£(x,y) 



where 



Dpi* tf m)v = m 2 <tf m)v , 



and the eigenfunctions h^(y) are determined from the equation 

h<"(y) - «(y) + e^y\m\h<(y) + 5{y)A l ih<{y) = 0. (43) 

In the region y > y in (22) the term with the increasing exponent becomes dominant. The approxi- 
mate metric is 



rf S 2 = (fy 2 + e 2 W ^_ df 2 + Y 7 ^ j Vabdx a dx b J ( J J ) 

The equation for the eigenmodes is 

- V£(v) + e^'m* /£ (y) = (45) 

Eqs. (43) and (45) are solved in appendix B. We show that the norm of the function is smaller 
than that h^. Effectively, this allows to neglect the contribution from the region y > y. Solving the 
equations we neglect t dependence of y(t) which in the final expressions can be effectively sent to 
infinity. In the radiation-dominated period of late cosmology e 2 ^ < 1 and (dy / dt) / fxy <C I. 
We obtain the spectrum 



m n ~ fie m 



(nn + |) (46) 
and the normalized eigenmode h m (Q) 

h m (0) ~ (/ie-^) 1 / 2 . (47) 



For the following we need the sum J2 n ^m„(0)> where m n is determined by (46). Because of a narrow 
spacing between the levels, we change summation to integration and obtain 

£C(0)^/^e-^/^. (48) 

The integral (48) is independent of y. The same measure of integration was obtained in [11], where 
the authors used the graviton modes of the RS2 model without matter, in which case the integration 
over y extends to infinity and the spectrum is continuous. Similarity of the results can be traced 
to the fact that we performed calculations in the period of late cosmology neglecting the terms of 
order 0(p//i) as compared to unity. This result can be anticipated also from the equation for the 
eigenmodes which does not contain matter terms. Time (temperature) dependence of the Universe 
enters only through the y , and calculations of appendix B are effectively performed by taking the 
limit y — > oo. 



6 Production of Kaluza-Klein gravitons 

Let us calculate the rate of production of Kaluza-Klein gravitons in interactions of particles in hot 
matter on the brane in the radiation-dominated period. 

Production of Kaluza-Klein gravitons in reactions of particles localized on the brane is calculated 
with the interaction Lagrangian 

I = kJ d 4 x^h^(0,x)T^(x), 

where T^ v is the energy-momentum of particles on the brane. The Boltzmann equation which 
determines evolution of the energy density pit) of matter on the brane is [20, 10, 12] 7 

f + ^=-EE/g/ ( 2 gU ^ WWM^ + *>-!>)■ (49) 

Here /' are the Bose/Fermi distributions of colliding particles and M l n is the amplitude of reaction 
+ ip l — > G, where and ip are standard model particles on the brane (vector, spinor, scalar) and G 
is a state of mass m n from the graviton Kaluza-Klein tower 8 . The squared amplitude of annihilation 
into graviton is [10] 

IKI 2 = A~h 2 n (0)s\ (50) 

where Ai = A s , A v , Af = 2/3, 4, 1 for scalars, vectors and fermions and s 2 = (k\ + k^) 2 . The sum 
over the graviton states is transformed to the integral following (48). After integration over the 
angular variables of the momenta of interacting particles Boltzmann equation takes the form 

d l + m (51) 



;(2vr) 4 J n 



/4(0)m 4 / dhdhfi^fihXh + ^O 1- 



1 



2hk 2 



7 In this section we use the unnormalized energy density. 

8 We have written the rhs of (49) without the factor 1/2 (cf. [12]), becauce in annihilation reactions ip l and 'tp 1 are 
different particles [20, 21]. 



Integrating over to, we obtain 

(2feifc 2 ) 1/2 



+ 4ffp=- ^|y J drum 4 j dhdfafikJfifaKh + fa) (52) 

o 

/t 2 A315C(9/2)C(7/2)T 8 



7T 3 2 8 

where the sum extends over relativistic degrees of freedom 

a = J2 A = it + 9f{1 - 2 ^ 7/2)(1 " 2_9/2) + 4 ^ (53) 

i 

In the high-energy period, when all the standard model degrees of freedom are relativistic, A = 166, 2. 
Eq.(52) has the same form as Eq.(18) and the expression in the rhs of (52) is the component T 4 of 
the bulk energy- momentum tensor [12] (see (69) in the next section). The physical content of both 
equations is that the energy loss on the brane is due to graviton production. 



7 The energy- momentum tensor in the bulk 

Corrections to the leading-order solution to the Einstein equation are expressed through the com- 
ponents of the energy-momentum tensor of the emitted gravitons in the bulk. In this section we 
determine the y-dependence of the components of the energy-momentum tensor in the bulk in the 
period of late cosmology. 

After being produced in reactions on the brane, gravitons move freely following the null geodesies 
in the bulk. Coordinates of particles on the geodesies are parametrized by an afhne parameter A. In 
the region < y < y we integrate geodesic equations using the approximate metric (33) 



Geodesic equations are 



ds 2 = dy 2 + e- 2 Wg^{x)dx?dx v . 



d 2 x a dx b dy dx b dx° _ 

~dv +2L b y^\d\ + 21 60 rfA 1\ ~ u (54) 

d ec „-^^n dec d%\ „.^^n dec dec „ , , 

Y +2Toy ^ o i +2Tab ^f=° (55) 

' / " // + r oof^V + r^^ = 0, (56) 



d\ 2 w \d\ J ab d\ d\ 
where connections are 

ya _ ya _ ra^ pO _ pO _ • y.y _ _ -2/ij/ y'y _ 2 . /p -2/ij/ 

1 by — °bh L J 1 bO ~ °b i L Oy ~ r"> L ab ~ <lab aa , 1 00 — P e ) 1 ab ~ U 'labH'V 

Integrating the geodesic equations, we obtain 

dx a C a e 2 ™ 



(57) 



dx = ^{* + i?fcr) (58) 



d\ a 2 (0,t) 

dt o 2 m ( n 2 , C a C a \ 1 



(|) 2 -C V^C, 



(59) 



in 



Using these relations, it is straightforward to verify that 

rlnr a d-r^ 

■ const (60) 



dt\ 2 dx a dx b 

Tx) + ar]ab lxlx 



is independent of A.. 

The energy-momentum tensor of the bulk gravitons is taken in the form of free radiation of 
massless particles [12] 

f ij (x) = J d 5 p^5(p t p l )f(x,pW, (61) 

where f(x,p) is the phase space density of the distribution function of gravitons. The function f(x,p) 
is determined by identifying the interaction term in the rhs of the Boltzmann equation (52) and 2T04 
in the {04} component of the Einstein equations (18) written in the radiation-dominated period, 
p = p/3, in terms of the unnormalized energy density p [12] 

^ + 4Hp = 2f 04 . 
dt 

The components of momentum p % of gravitons subject to condition pip 1 = are proportional to 
the tangent vectors to null geodesies 

(n na \ V 2 fia iiy 

C ° 2 + 'aWT)J ' P 4 = Ar (^K^, P a ^ = N (y^^Ty ( 62 ) 

where N(y, t) is a normalization factor determined below. Integrating the geodesic equations (57)- 
(59) with the initial conditions 

y (A = 0)=0, t(A = 0) = * , x a (\ = 0)=x a , 
we obtain x l = x' l (\,C\t ,x a ). Integrating (59), we have 



f 'dy'ic + cy^'r'/^x. 

Jo 



From this relation is determined A = \(y,C % ). Inverting this relation, we obtain y = y(\,C l ). 
Because of the initial conditions, y(0,C l ,x a ) = 0. 

The phase space density of non-interacting massless particles f(x\ p l ) as a function of A is constant 
along a null geodesic 

d\ ~ P dx* dp kKsPP " U - [M) 

Because the function f(x l (\, C l , t , x a ), p l (\, C\ t , x a )) = /(A, C\ t , x a ) is constant along geodesies, 
we have 

f(x\\, C\ t , x a ), p\\ C\ t , x a )) = /(0, t , x a , p l (0, C\ t , x a )) = f(C\ t , x a ). 
Let us consider the {04} component of the energy-momentum tensor 

f°\y,t,x a ) = J d 5 p^5( Pl p l )f(x,p) P y. 

1 1 



Changing integration variables p % to C\ we obtain 

d(p°,p 4 ,p a ) 



f 04 (y,t) = J d 5 C l e-^ y a 3 (0,t) 



d(C°, C 4 , C a ) 
N 2 (y, ty™ [c, 2 + (C 2 e 2 ^ + C 4 ) 1/2 . 



5(N 2 C 4 )f(C\t ,x a ) (64) 

1/2 



The Jacobian of the transformation is 



/iVe 2 ^ 3 Ne 2f * y C N 

U(0,0J (Co 2 + C a CVa 2 (0, t)) 1/2 2 (C 2 e^y + C 4 ) 1/2 ' 



(65) 



We traced the distribution function along the geodesies to A = 0, i.e. to the hypersurface y — 0. 
Separating the ^/-dependence, we have 



T %,t)=T> ; t)(^|ye 6 -. (66) 

Following [12], we expand the graviton momentum p % in the orthonormal basis adapted to the foliation 
of the space-time 

u* = (1, 0, 0, 0, 0)e w , n* = (0, 0, 0, 0, 1), e» = (0, 1 , 0, 0, 0)^y. (67) 

We have 

p' = Eu { + mrC + q a e i a , 

and pip 1 = —E 2 + m 2 + q a q a . 

Projections of the components of the bulk energy-momentum tensor on the plane y = 

T M = TVnj, f m = f ij u iUj , f u = f ij n inj (68) 

have the form [12] 9 

where T is the temperature of the Universe at the time t. 
Separating in Tgi the y- dependence, we have 



(70) 



The same relations hold for all components: fjj(y,t) = (N(y, t)e m /N(0, t)) 5 T-j(0, t). 

Our next aim is to determine the normalization factor N(y, t). Lacking an independent definition 
of momentum of massless particle which could enable us to fix the normalization factor, let us consider 

9 Cf. (52). 
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the energy-momentum tensor conservation equations in the bulk V {T 1 ^ (y , t) = and try to extract 
from these equations an information on y-dependence of the energy-momentum tensor. 
The component j — of the conservation equations is 

ViT^y, t) = «9 4 T 40 + (3r° 4 + H 4 )T 40 + d f 00 + (2r° + T a a0 )f 00 + T° ab f ab = 0. (71) 

Substituting the explicit expressions for the Christoffel symbols, we obtain 

;i / . _ j / 11 :' ■> " i o a \ rpAO , ;. •/•Ill) . / ., " . .y"\ will) . "" .. rpab 



ViT°(y, t) = <9 4 T 40 + ( 3- + 3- ) T 4U + <9 T UU + ( 2- + 3 7 ) T uu + — Vab T ab = 0, (72) 

where 



n a I \ n a / n 2 



T 04 e-™ = To4, f 00 e-^ y = f m , f 44 = f u , T ab e- 2 ^a" 2 (0, t) = f- al . (73) 

In the region < y < y, using the approximate metric (25), we have 

n'/n ~ — /i, a'/a ~ — /x, a/a~H(t), h/n ~ H(t)p/ p. 

Taking into account that in the radiation-dominated period T ~ a(0, t)" 1 , we obtain that <9 Too(0, t) ~ 
-8H(t)f- - (0,t). 10 Thus, Eq.(72) can be written as 

(9^ - 6/i^)T 4(0,t) + i/ e ^(-5Too(0,t) + r]a b T- ab -(0,t)) = 0, (74) 

where tp = e^ y (e m N(y,t)/N(0,t)) 5 . The trace of the energy-momentum tensor of massless particles 
is zero, = 0. Using (68), we obtain a relation 

TabVab = TOO ~~ ^44 ■ 

Eq.(74) can be rewritten as 

d y ip - 6pip + ipHe^B = 0, (75) 
where B = (AT m + f u )/(-f u ) ~ 6.44. Solving (75), we obtain 

Next, let us consider the j = 4 component of the conservation equation of the energy-momentum 
tensor 

v,r 4 ( y , *) = d A f 44 + (r° 4 + r^ 4 )T 44 + d f 04 + (r° + r a a0 )f 04 + r 4 T 00 + r 4 ab f ab 



d A f 44 + ( - + 3- ) f 44 + d f 04 + ( - + 3- ) f 04 + n'nf 00 - a'a^f ab = (77) 
n a J \n a 



dy<p - hpup - <pHef iy = 0, (78) 



Repeating the same steps as for the component j = 0, we obtain 

, 57-54(0, t) 
f u (0,t) 

where ip = (N(y, t)e^/N(0, t)) 5 and B' = -5T M /f u ~ 6.21. We have used that at y = there is a 
relation 

n 'nf 00 - aai] ab f ab = -p(f m - i] ab f ab ) = -pf U - 

Solving (77), we have 



N(y,t)Y f B'H(t) 



= exp <^ ^(e™ - 1) } (79) 



N(0,t)J * I /x 

It is remarkable that the constants B and B' which appear in different equations are close to each 
other. The small difference between the constants can be attributed to approximations to the exact 
equations. 



10 We neglected a weak T dependence of g*{T). 
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8 Einstein equations with the bulk energy-momentum ten- 
sor 

Having the expressions for the bulk energy-momentum tensor T l \ we are in position to find the range 
of t and y within which solutions of the Einstein equations with the bulk energy-momentum tensor 
included are sufficiently close to the leading-order solutions, i.e. can be considered as perturbations. 
In particular, in Sect. 4 we determined the spectrum of Kaluza-Klein tower of gravitons in the leading 
approximation which is possible, if p/p <C 1 and if corrections to the leading approximation are small. 
Let us consider the ([]) component of the Einstein equations in the bulk 

a"(y,t) a' 2 (y,t) a 2 (y,t) n 9 2k, 2 ~ n . . , . 

a{y,t) a 2 (y,t) n 2 {y,t)a 2 (y,t) 3 

The ratio of the bulk term to 2p 2 is 

K 2 f$(y,t) _ K 2 f m (y,t) _ ( N(y,t) \ 5 n 2 f m (0,t) 



_ 2 ^4 180p 



(y) e WV ^ (SI) 

3/i 2 3/i 2 \N(0,t)J {V) 3p 2 1 } 

where 

K 2 f m (0,t) 21AC(9/2)C(7/2) / 180 \ 2 (p V 2 _ /p V 2 
Here we have used (19) to write 

where g* — g s + g v + 7g//8. The coefficient varies in the range from K ~ 0.2 at characteristic 
energies of the nucleosynthesis, T ~ 1(H 3 ^ 4 ) GeV, to K = 0.025 at the energies ~ 10 2 : 3 GeV at 
which all the standard model degrees of freedom are relativistic. In this range the ratio p/p changes 
from icr 22 ^ 29 to lO^ -3 ^ 1 ^. Expression (81) can be written as 

^(y,t) =K (A\s„lHB {e „_ l) \ 



3p 2 \p ) \ p 

where HB/p ~ 9.1(p/ p) 1 ^ 2 . In the region y < y, where e 2 ^ = 1 + 2p/p, for p/p <C 1, we have 



K 2 Tg(y,t) <CK [P 



2 



3/i 2 \p 

where C = 0(1). Fixing ti so that p(ti)/p <C 1, the ratio (81) can be made arbitrary small. 
Integrating Eq.(7), we obtain the (04) component of the Einstein equations as 

a(0,t)exp - / dy . (84) 



»(//•') ' V ./o 3a 

Using the formulas of Sect. 7, we have 



f 

Jo 3a ' J Ml 

K 2 ni(o,t) r w r ^ 



3/7 



jT e w exp |-^(e"»' - 1)| dy'. 
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In the period of late cosmology, using the approximate metric (25), we can set a(y', t)/a(y', t) ~ H(t). 
The factor at the integral (85) is 



iftUM 3i5A a fp(ty\ 3/2 = fp\ 3/2 

At the temperatures of the Universe of order 10~ 3 GeV, we have A = 18, g* = 10.75 and K ~ 0.12. 
At the temperatures at which all the Standard model degrees of freedom are relativistic, K ~ 0.011. 
The integral in (85) is 

/ = y dssV 6 *- 1 ' = e 6 £(/(fc) - /(fce^)), 

where b = HB/p and 

f(b) = e~\l + b + b 2 /2\ + b 3 /3\). 
For 6 < 1 approximately /(&) ~ 1 — b 4 /12, and for y < y the integral is estimated as 



/4(e--l)<2^exp{ 9 ,(£) 1/2 } 



2 

For small enough p(ti)/p the expression in the exponent in (84) 



*u> ^ru 1 > (87) 

can be made arbitrary small. 

To conclude, contribution to the energy-momentum tensor due to graviton emission modifies 
solutions of the Einstein equations, however, in the period of late cosmology corrections to the 
leading-order solutions are numerically small as compared to the leading-order terms. 

Let us consider the Friedmann equation (16) with the terms %j included. Calculating the first 
integral in (16), we obtain 



2k 2 

'04 



3a 4 (0,t) 



/ dt'f Q4 {0,t'){^ + p{t'))a-\0,t') = (88) 
Jti 



12 \fiti fit J 288 \(pti) 3 (A**) 



where we have substituted (69) 



2 



«V M = W = _ 315.4C(9/2)C(T/2) M> ^ 



2 9 7T 3 \g*lT 2 

The integral (88) has a strong dependence on the value of the lower limit. Therefore, in the integrand 
the slowly varying functions can be taken at the times when all the Standard model degrees of 
freedom are relativistic. In this period A = 166.2, g*{T) = 106.7 and A^ ~ —0.126. In the 
leading approximation p(t) ~ l/8pt 2 and 1/ pU ~ {8p{ti) / p) 1 / 2 Taking p{ti)j p ~ 0.1 4- 0.001 and 
p ~ 10- 12 GeV, we have 1/pU ~ 0.9 -=- 0.09 and T x ~ (5.1 -=- 1.6) • 10 2 GeV '. For p(t t )/p ~ 0.1 we 
obtain 

7 04 ^ -2/ip • 0.0048. (90) 
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In the same way, calculating 7 44 , we obtain 



'44 



3<" 4 (0,«)7,, 
where 

^ = . W = 3^(9/2)C(7/2) / 180 V 2 ^ g jQ 3 

47T 4 \g*ir 2 J 

This yields 

7 44 ~ 2/ip • 0.00094. (92) 

9 Graviton emission in the bulk and nucleosynthesis 
9.1 Graviton emission in the period of late cosmology 

Let us compare abundances of 4 He produced in primordial nucleosynthesis calculated in the models 
with and without account of the graviton emission in the bulk. We consider the radiation-dominated 
period of late cosmology in which production of gravitons is sufficiently intensive. In the period of 
late cosmology the leading approximation of the Friedmann equation coincides with the standard 
cosmological model. 

We estimate the effect of the graviton emission on nucleosynthesis by solving Friedmann equation 
perturbatively keeping in the generalized Friedmann equation (16) the term linear in the radiation 
energy density and the terms containing the bulk energy-momentum tensor. The Weyl radiation 
term and the term quadratic in radiation energy density can be treated perturbatively also. We have 

H 2 ~ 2/ip - J 04 - J44, (93) 

where the integrals /04 and 7 44 were defined in the previous section. We use also the 5D energy 
conservation equation (18) which is written as 

p + AHp = -^f. (94) 

Let p and H be the energy density and the Hubble function calculated in the model without inclusion 
in the Einstein equations the terms with the bulk energy-momentum tensor and the Weyl radiation. 
Defining 

p = p + p 1 , H = H + H U 

where H 2 = 2pp and p\ and H x are perturbations, and separating in (93) and (94) the leading-order 
terms, we obtain 

2HH, = 2p Pl - y y dt' [A 0A (p + p) + A U H] p(t') (95) 
Pl +AHp 1 +ApH 1 = ^^. (96) 

Substituting in (95) p(t) ~ l/8pt 2 and performing integration, we have 

HA U (\ 1 



it u\ V HA GA 



11 1/11 

- - + 



ti t 24/i 2 \tf t 3 



192/i 2 \tf t 2 



(97) 



16 



Substituting expression (97) for Hi in (96) and noting that H — l/2t, we obtain 



Pi + -Pi 



A 



04 



192/z 2 t 3 

Solving this equation, we find 
A 04 



1 1 



Pl(0 = ^ + 



t 3 192/i 2 t 3 



— + 



24/z 2 



t-U 

t 3 



1 



1 



+ 



A 



44 



768/i 3 t 3 V*; 



1 

t 2 



1 



+ 



44 



768/x 3 t 3 \tj t t + t 



(98) 



The constant C\ should be determined by sewing solutions of Friedmann equation in the periods of 
early and late cosmologies. For a moment we set C\ = 0, i.e. look for a contribution from the period 
of late cosmology. For Hi we obtain 



Hi(t) 



A 



04 



1536/x 3 t 2 \t 



1 

T 2 



A 



44 



192/i 2 



1 1 



A 



04 



A 



41 



-V 



192 fiHf 



(99) 



Expressions (98) and (99) show that corrections to the leading terms are small, i.e. the perturbative 
approach is justified. Note that the leading term proportional to A^/ fJ>U was canceled in Hi. 
The mass fraction of A He produced in primordial nucleosynthesis is [20, 21] 



X 4 



2(n/p) 
(n/p) + 1 



where the ratio n/p is taken at the end of nucleosunthesis. Characteristic temperature of the onset 
of the period of nucleosynthesis (freezing temperature T n of the reaction n <->■ p) is estimated as the 
temperature at which the reaction rate ~ GpT 5 is approximately equal to the Hubble parameter [20] 

G F T n 5 ~ H. 

The difference of the freezing temperatures in the models with and without the account of the 
graviton emission is 

ST n H 1 lp(t n ) ( A. 



T 



5H 



04 



+ 



A 



44 



where 



P(tr 



2/i v l9 V*f 24 pt 



ln 3 g*(T n ) T* 



(100) 



(101) 



H 180 /iM 3 ' 

Substituting T n ~ l^GeV and /iM 3 ~ {pM pl ) 2 ~ 10 u GeV 4 , we find that the ratio 5T n /T n is very 
small. The equilibrium value of the n — p ratio at the freezing temperature 



cxp 



T 

± n 



is very sensitive to the value of T n . Substituting 5{n/p) n = (n/p) \n(n/p) n ST n /T n , we obtain variation 
of X 4 under variation of the freezing temperature 



5X A 



(n/p + 1) 



In 



6Tn 

T 



(102) 



which is also a very small number. 
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9.2 Graviton emission in the period of early cosmology 

Next, we make an estimate of the variation of X A due to the graviton emission in the bulk in the 
period of early cosmology, in which the p 2 term in the Friedmann equation is dominant. For the 
value of /i ~ 10~ 12 GeV, the characteristic temperatures of this period are above 5 • 10 2 GeV. 

To make an estimate of the graviton emission, we assume that the collision integral in the Boltz- 
mann equation and the expressions for obtained for the period of late cosmology remain qualita- 
tively valid in the early cosmological period. The new phenomenon in the early cosmological period 
is that some of the emitted gravitons can return to the brane and be again reflected in the bulk with 
a different momentum. These gravitons do not contribute to the component T 0A , because they are 
not produced, but reflected, but contribute to the component T44. The Friedmann equation contains 
a new term f$ representing the energy-momentum tensor of gravitons bouncing back to the brane 
[10,12] 



H 2 {t) = p 2 {t) + 2pp{t)+ pp w {t) 



(103) 



2k 2 



3a 4 (0,t) 



dt' 



f 04 (t')(p(t') + p) + f M (t')H(t') - TS(t')H(t') a 4 (0, t' 



Here the initial time t c is the time of reheating. Numerical estimates [12] and considerations from the 
Vaidya model [11] suggest at t ^> t c the dominant contributions from the integrals of the components 
Tij and the integral of the term fj$ are mutually canceled. 

Let again p and H be the energy density of matter on the brane and the Hubble function calculated 
in the model without the graviton emission in the period of early cosmology. We define 



P = P + P2, 



H — H + Ho 



where in the period of early cosmology p(t) ~ l/4t, and H ~ p{t) (see (27). It is assumed that p 2 
and H 2 are perturbations to the leading terms. From the Friedmann equation we have 



2HH 2 ~ 2pp 2 dt' \A 0A {p + p) + A 44 H - K 2 fSH p{t'). 



(104) 



The conservation law yields 



p-2 + 4:Hp 2 + 4pH 2 



A 04 p 



7$ 



(105) 



In the same way as in the period of late cosmology, we obtain the equation for p 2 

A AA ( 1 1 



2 1 
p 2 + -p 2 -- 



4 , t 1/1 1 

u In h 

12 I 1 



t c 4 V t c t 



+ 



48 \t c t 



(106) 



where we substituted p ~ l/4t. Integrating (106) with the initial condition p 2 (tc) = 0, we obtain 

' 2 ^ /I- 1 /i / Aaa f\ tr 



( ) _ A m p f t 1 t 2 c \ A 04 + A u t c 
p2{t) -^TV n Z~2 + 2^ + 9Qt r l 1 "^/ 48 W t 2 



(107) 



-^f* dyy£ dxKtffflix) 
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The time t c of reheating is estimated for the reheating temperature Tr ~ 5 ■ 10 6 Gel / [23]. From the 
relation p(t c )/p ~ \/Apt c we have 

1 ^ l67r 3 g*(T R )Tx . 



/it c 90M 3 



~ 4 • HT, 



where we have put M ~ (/iM 2 ,) 1 / 3 ~ 10 8 Ge\/ and g*(T R ) ~ 10 2 . It follows that i,/t c ~ l//it c ~ 10 16 . 
In (107) there is a large term (A 04 + A 44 )/i c . Omitting the small terms, we have 

IK) 24 V ^ 2 J m c 48 t 12t 2 J tc yy J tc 44 KJ 1 ; 

On dimensional grounds the term K 2 T^(t) at small t has the following structure 

^fi?(t) = | + ^ + ^ + ... (109) 

Performing integration of the last term in (108) and taking t ~ ti, we obtain 

/ , \ A 04 p ( ti 1\ A AA A 04 + A M b 2 b 2 hp ti 

poiti) ~ In In— 110 

l2KU 24 V tc 2 J 48t, 96t c 24t c 12*, 24 t c K J 

Next, we equate P2{U) and Ci/tf in (98). At the times t ~ U : where pti ~ 1, we have 

^ V 24 V Pt c 2 J 48 y 96t c 24t c 12 24 /it c V ; 



The term d/t 3 is 



Ci 1 
t 3 ~ /x 2 t 3 



, | ^1 /'in J_ 1^ _ Au ^\ _j_ ^04 + ^44 _ _&2_ &1 ^! n 1 

; ' 1 24 \%t c 2/ 48/ 96t c 24t c 12 24 H pt c 



(112) 



In the period of late cosmology the term C\/t 3 generates in Hi an additional contribution AHi 
pApi/H 



' ' 1 24 \%t c 2 ; 48 7 96t c 24t c 12 24 %i c 



AHi = lQp 

where p(t) ~ 1/ 8yut 2 is the energy density in the period of late cosmology 



(113) 



If the term (A 04 + A 44 )/96t c ~ 3-10 A /t c was not canceled, it would produce in pi the contribution 

A m + A 44 1 3 • 10" 4 



Ap 

From (97) we would have 



1 " 96t c (/it c ) 3 (pt)H c 



At time of the nucleosynthesis 



AH i A* a 1.2 -10- 3 

1 4tt 3 ^(T to ) g 

8(/xt n ) 2 90 KO 2 ' 
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For T n ~ 10~ 3 GeV, we have pt n ~ 10 12 . From (114) we obtain AH^tn)/ H(t n ) ~ 4, which is too 
large a value, and would contradict the experimental data. 

Assuming that the large terms in (113) are canceled, we have 



Because at the period of nucleosynthesis AHi/H ~ (p(t n ) / p) 1 ^ 2 is a small number, contribution from 
the early cosmology would result in a small variation of 5X 4 /X 4 . 

10 Conclusions and discussion 

Calculations of this paper where performed under the restriction to the period of late cosmology, 
when p(T)/p < 1, where p(T) is the normalized radiation energy density of matter on the brane, p 
is the scale of the warping factor in the metric. For p ~ !0~ 12 GeV, which we adopted in this paper, 
the limiting temperatures of the Universe at which the approximation of late cosmology is valid are 
of order T) ~ 5 • 10 2 GeV . 

In the period of late cosmology it was possible to make a number of approximations, which 
enabled us to obtain the analytic expression for the energy loss from the brane to the bulk due to 
graviton emission. We calculated the spectrum of gravitons and found the energy-momentum tensor 
of emitted gravitons in the bulk. The Einstein equations were solved perturbatively, taking as the 
leading-order solution that without the graviton emission. Using the expressions for the energy- 
momentum tensor of gravitons in the bulk, we estimated corrections to the leading-order solution to 
the Einstein equations and showed that the perturbative approach is justified. We estimated also 
corrections to the leading-order terms in the generalized Friedmann equation and showed that in the 
late cosmological period they are sufficiently small as compared to the leading-order terms. 

Graviton emission changes the cosmological evolution of matter on the brane and thus production 
of light elements in primordial nucleosynthesis. Solving the system of the generalized Friedmann and 
the 5D energy conservation equations, which included the terms representing the energy flux from 
the brane to the bulk, we found the difference of the mass fractions of A He produced in primordial 
nucleosynthesis calculated in the models with and without the graviton emission 



Here p is is the radiation energy density on the brane in the leading approximation, without account 
of graviton emission, T n ~ 10~ 3 GeV is the freezing temperature of the reaction p -h- n, and a < 1 
is a small number. The ratio p(T n )/p ~ 10~ 26 is a very small number, and thus the difference of 
abundances of 4 He calculated in both models is small. Crude estimate of the contribution to 6X4/ X± 
from the early cosmological period based on the assumption that in the Friedmann equation the 
large contributions from the terms containing the energy-momentum tensor of the emitted gravitons 
cancel due to the bounce of the emitted gravitons back to the brane [12] again shows that it is a 
small number. 

In papers [10, 12] as a measure of energy density loss on the brane due to graviton production 
was introduced the integral 




(115) 





ti 



where — Ap is the collision integral in the rhs of the Boltzmann Eq.(52) which gives the rate of 
graviton emission. Because in [10, 12] and in the present paper the collision integrals are identical 
(up to the factor 1/2), expressions for Qi ost up to the factor 1/2 are the same. Substituting expressions 
for Ap and p, we find 

a 04 n r 



lost 



24/i \ti t 

where t\ is the limiting time at which approximation of the late cosmology is valid, pi\ ~ 1, 
and A 04 = —0.126. However, calculating correction Hi(t) to the leading-order Hubble function 
H(t) = l/2t, we found that in Hi(t) the terms H(t)A 04 /pti are canceled. The ratio of the would-be 
correction H(t)A 04 /ASpti to the actual correction p(t)(A 04 / '192 p?t 2 + A 44 /2Apti) is of order (p/p) 1//2 . 
As discussed above, at the temperatures of nucleosynthesis this is a small number. 

An alternative approach to calculation of emission of the Kaluza-Klein gravitons from the brane 
to the bulk was developed in [13, 14] and refs. therein. In this scenario there are a visible brane 
moving in the AdS bulk and a static brane. The motion of the visible brane is determined by 
evolution of matter on the brane. The moving brane acts as a time-dependent boundary to the 5D 
bulk leading to production of gravitons via the Casimir effect. Although this picture in appearence is 
very different from that of the present work, many of the intermediate formulas and final expressions 
are rather similar. In particular, the energy density of the emitted Kaluza-Klein gravitons on the 
brane was calculated to be pkk ~ l/° 6 ? where a(t) is the scale factor in the induced metric on the 
brane. The energy density loss on the brane is dpKx/dt ~ a/a 7 . In the radiation-dominated period 
of late cosmology, using the Friedmann eqation (a/a) 2 ~ 2pp, it is obtained that T ~ const /a(t) 
[20, 21] and dpxK/dt ~ T 8 in correspondence with the present work and [10, 11, 12]. 

Graviton emission in the bulk was estimated also in the framework of DGP-like models. In the 
model with the action [24] 

S = M 3 j d 4 x J* dyy/-g&R® + M 3 r c j d 4 x^-g^R^ 

where M is in the TeV range, M 2 t = M 3 (R + r c ), the ratio of the rate of the graviton emission in 
the bulk to the cooling rate of matter due to cosmological expansion is 



dp 
~dt 



dp 
~dl 



T 

-a, 



exp M p . 



where a ~ 10 4 . 

In [25] was proposed a model with the action 



S = Ml J d A xdy^g~^R^ + M 2 J d A x^g~^R^ + J d 4 Xv ^g~^L S M, 

where M* ~ l0~ 12 GeV. The strong gravitational coupling, 1/M 2 , presumably, is renormalized by 
the loops of the standard model fields down to 1/M 2 t The ratio of the rate of the graviton emission 
to the cooling rate of matter was estimated as 



dp j dp 
dt J dt 

This ratio is less than unity for T < 10 2 GeV. 



rp3 

exp ~ M~M' 
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A 3D brane in 5D bulk 

From the system of 5D Einstein equations and junction conditions follows the 4D form of the Einstein 
equations 

(^(4) _ 2 T (4) 

where the Einstein tensor G$ is formed with the 4D metric g^, and Explicitly [5] 

k\T$ = pK% u + + ^- (f^ + (f u - -T^j - E^, (116) 

where we have set a — p. Here is quadratic in r MI/ , and is the 'electric' part of the 5D Weyl 
tensor [5, 7]. The component T^q is 

n\T^l = -6pp - 3p 2 - pp D: (117) 

where 

2k 2 /-„ . T\ n 

Noting that = — 3iJ 2 , we obtain the generalized Friedmann equation 

if 2 = 2^ + p 2 + ^. (118) 

From the 4D Bianchi identity D^G {4) ^ = follows the 4D conservation law D^ 4 ^ = 0. The 
zero component of the conservation law, D^T^^ = 0, is 

T (4) ° + 4#T (4) ° - #T (4) £ = 0. (119) 

Substituting in this equation T^q from (117) and expressing p via (18), we obtain (cf. [9, 12]) 

2k 2 

Pd + AHp D + [T 04 (/i + p) + Hf u ] = 0. (120) 

p 

B Eigenmodes and eigenvalues 

Solution of the equation 

hl{y) - Ap 2 h m {y) + e 2 ^m 2 h m (y) = (121) 

is 

hm(y) = CMme^) + C 2 N 2 (me^), (122) 

where 

m 

m = — . 
p 

The term with 8(y) is taken into account by the boundary condition 

'dhjjj) 2 ] 

— j h —h m {y) = 

L m J y=0+ 

which yields the relation 



Using this relation we obtain the eigenfunctions in the regon y < y the form 

h^y) = C [(JVi(m)J 2 K l1 ' 1 ) - Ji(m)N 2 (me M )] . (123) 
11 Introducing the functions 

f k (y) = N 1 (m)J k (me' i M) - J^N^me^) 

and using the formula 

f x 2 

/ dxxZl(ax) = — \Z 2 n {ax) - Z n+1 (ax)Z n _ 1 (ax)], (124) 

where Z is any linear combination of the Bessel functions, we obtain 

y 

ll^ll 2 = 2C 2 / e^fi(y)dy = — [e 2 ^ (f 2 (y) - h{y)fM) - {fW) - /i(0)/ 3 (0))] . (125) 
o 

Typical masses (energies) of the emitted Kaluza-Klein gravitons are of order of temperature of the 
Universe T. In the case of fi ~ 10~ 12 GeV we have m/fi ^> 1. 

For large values of arguments of the Bessel functions, we obtain 

fi(y) ~ ( — —J [— Ni(m) cosz + Ji(m) sinz] (126) 
/ e-wx 1 / 2 

f 2 {y) — ( — — I [NAm) sinz + JAih) cosz] 
\ irfh J 

/ e -/iJ/\V2 

/3 (y) — ( \ [Ni (to) cos z - Ji (m) sin 2;] , 

where z = roe™ — tt/4. For large to, substituting asymptotics of the Bessel functions, we have 
J 2 (rh) + N 2 (rh) ~ 2/nm. At the upper limit of integration we obtain 

f2(y)-fMh(y)^ (4-) e-'*. 

At the lower limit we have 

/ 2 2 (°)= f4V, /i(0)/ 3 (0) = 0. 



Combining the above expressions, we have 

.<m*_cv_h_V 



ik" -tUJ (e " _1) ' (127) 

In the region y > y in (22) the term with the increasing exponent becomes dominant. We obtain 
the equations for the eigenmodes 

h^Uy) - ^ 2 h<(v) + e- 2 ^mlh> m> {y) = (128) 



L To simplify the formulas wefrequently omit the sub (superscripts) > and <. 
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with solutions 

K>(V) = CMrhe-w) - C 2 N 2 {me-™). (129) 

For large y, such that me _ra <C 1, the function iV 2 (me _M?/ ) ~ (me - ' 1 ')" 2 rapidly increases, and to 
have normalizable eigenfunctions we set C 2 = 0. The norm of the function h^(y) is 



e- 2m J%(me- m )dy = C 2 ——[J 2 (me~^) - J^me"" 5 ) J 3 (me-™)]. (130) 

y 

At temperatures T >> /i at which production of gravitons is sufficiently intensive me"™ ~ (T/ /x)(p/ /i) 1//2 
1, and we can substitute the asymptotics of the Bessel functions 12 . 
Asymptotic of the eigenfunction (123) is 

1C 

h>(y) = — ^e~™ /2 cos(m(e" s - 1)). (131) 
Instead of sewing the oscillating functions h m (y) and h^(y), we sew the envelopes of their asymptotics 

^ e -^ = a,(^y /2 (132) 

7im \ Tim J 

giving 

/ 2 \ 1/2 - 
C 1 = C[ e~ m 



\irrh J 



Using this relation, we obtain 



Because > 1, the norm (133) is smaller than (127). 

Effectively, we neglect the contribution from the region y > y and impose the condition h m (y) = 0. 
Substituting h m (y) from (131), we have 

cos(m(e MS - 1)) = 0. 
From this equation we obtain the spectrum 

m n ~ lie'™ (nir + ^ . (134) 

Normalization constant C determined from condition \ \h m \ \ = 1 with the norm (127) is 



C ~ ^^V2 e -w/2. (135) 



The normalized eigenmode /i m (0) is 



0/7 

MO) = -4 ~ (pe-^) 1 / 2 . (136) 
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12 Taking T ~ W 2 GeV, we obtain from (19) p/fi ~ 1(T 4 and T/fj, ~ 10 10 . 
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